Thermodynamics of QCD low-energy models 
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We study the equation of state of a strongly interacting theory of relativistic bosons and chiral 
fermions in the vicinity and above the chiral phase transition temperature. Our model resembles 
presently used low-energy models of QCD in many ways, but its simplicity allows us to study 
systematically various approximation schemes by means of a derivative expansion. In particular, we 
compare the results for the phase transition temperature, equation of state and the thermal mass 
of the scalar fields obtained in various approximations. We find that the large- N c approximation, 
being the zeroth order of our approximation, deviates significantly from approximation schemes in 
which effects of (Pseudo-) Goldstone modes have been taken into account, even at high temperatures 
and even if we allow for a finite explicit symmetry breaking. The important role of the (Pseudo-) 
Goldstone modes is also manifested in the phenomenon of "local order" which is missing in large- N c 
studies. 

PACS numbers: 11.10. Wx,11.30.Rd,12.38.Aw,64.60.ae 
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I. INTRODUCTION 

The thermodynamics of strongly interacting relativis- 
tic bosons and fermions play a prominent role in various 
fields of physics, such as heavy-ion collision experiments. 
There the theoretically computed equation of state of 
QCD is used as an input in order to describe the expan- 
sion of the so-called fireball [1]. However, the computa- 
tion of the equation of state is a highly- involved problem 
since the high-temperature limit is not a weak-coupling 
limit, at least for the soft modes, see e. g. @, Q. On the 
contrary, by increasing the temperature we interpolate 
continuously between the 4d theory and the underlying 
3d theory which does not generally have to be weakly 
coupled. Indeed, QCD in 3c? as well as simple scalar field 
theories in 3d have shown to be strongly coupled. For 
this reason, a naive perturbative expansion of the parti- 
tion function at high temperatures is bound to fail. 

Several methods have been developed to study the 
equation of state. Lattice simulations are definitely the 
most powerful tool that we have currently at hand for 
a computation of the equation of state of strongly cou- 
pled theories. However, since lattice simulations are per- 
formed in a finite volume and at finite explicit symmetry 
breaking (i. e. current quark mass in the QCD con- 
text or external magnetic in the context of Ising models), 
the result for the equation of state has to be extrapo- 
lated to the continuum limit. The presence of dynami- 
cal fermions complicates lattice simulations further. On 
the other hand, continuum approaches have to deal with 
complementary problems, such as necessary truncations 
of the action of the theory. 

One prominent approach to a systematic study of 
the thermodynamics of strongly interacting theories 
are hard-thermal loop resummation techniques, see 
e. g. Refs. SSEHSiMEl- For exam P le , it has been 
shown in this framework that the momentum dependence 



of the n-point functions longitudinal and transversal to 
the heat-bath has to be resolved carefully in order to 
properly describe the equation of state of hot QCD above 
the phase transition. On the other hand, low-energy 
models of QCD based on the NJL model [12| have been 
extensively used to study the thermodynamics of strongly 
interacting hadronic matter in the vicinity of the phase 
transition. However, these models are not confining and 
do not include gluodynamics. Therefore phenomenolog- 
ical extensions of these models have been proposed, see 
e. g. [H 03, [H [H, 03 E 03 HI. By means of this, 
the thermodynamics at high temperatures are corrected 
in the sense that the missing gluon contributions are now 
incorporated. Still, the models are mostly studied in the 
limit of many colors - the large N c limit. As a conse- 
quence the contributions of the scalar fields, i. e. pions 
in QCD, are not properly taken into account. Neglect- 
ing the fluctuations of the scalar fields affects universal 
aspects of the theory in the vicinity of the phase tran- 
sition, such as critical exponents, but also non- universal 
properties, such as the phase transition temperature or 
the equation of state. 

In the present paper, our aim is to study systemati- 
cally the role of the neglected fluctuations of the scalar 
fields and of the non-trivial momentum dependence of the 
propagators of the bosonic and fermionic fields. To this 
end, we apply a derivative expansion to a strongly inter- 
acting theory of bosons and fermions with iV c colors but 
only one quark flavor. The low-temperature limit of our 
model is dominated by the presence of a single Goldstone 
boson (pion) rather than three as in QCD with two quark 
flavors. However, in the spirit of QCD low-energy models 
our model is as close as it can be to the large- ./V c approx- 
imation, where loops involving scalar fields are dropped. 
As we shall discuss below, the large-iV c approximation 
represents the lowest order of our derivative expansion of 
the effective action. 
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The paper is organized as follows: In Sect. |TT] we in- 
troduce our model and the non-perturbative framework, 
namely the functional Renormalization Group (RG) [HJ , 
which we employ to compute the thermodynamic prop- 
erties of our model. The results for the equation of state 
and the thermal mass of the scalar fields are then given 
in Sect. Mil In addition we discuss the dependence of the 
equation of state on explicit symmetry breaking (i. e. 
current quark mass) and on the ultraviolet cutoff. Our 
concluding remarks, including a discussion of possible fu- 
ture extensions, are presented in Sect. IIVI 



II. FUNCTIONAL RENORMALIZATION 
GROUP AND LOW-ENERGY MODELS 

A. The model and the studied truncations 

Throughout this paper we work in d = 4 dimensional 
Euclidean space and employ the following ansatz for the 
effective action for a strongly-interacting theory of bosons 
and fermions: 

T = J d 4 x <.Z^4> a i$ip a + -^=tp a (f ■ $)^a 

+ Iz (9 AI <f>) 2 + C/($ 2 )| , (1) 

where the index of the fermion fields denotes an inter- 
nal index, e. g. a color index in QCD, and runs from 
a = 1,...,N C . The scalar fields are combined in the 
0{2) vector <I> T = (a, ir) and we have used f — (75, i • 1) 
in order to define a chirally invariant Yukawa interac- 
tion. We assume that the bosons are composite degrees 
of fermions and do not carry an internal charge, e. g. 
color: a ~ ipa^pa and ir ~ ipa"f5ipa- In the spirit of 
hadronic models, our ansatz for the effective action rep- 
resents a low-energy model for QCD with one quark fla- 
vor. However, we rush to add that our interest is not the 
study of one-flavor QCD. In this case, we would have to 
take the effects of topologically non-trivial gauge config- 
urations into account which are prominent in one-flavor 
QCD since they induce mass-like fermion interactions 
and break the U A {1) symmetry [H SI SI SI . Here, 
we are interested in an investigation of the effect of light 
bosonic degrees of freedom on the thermodynamics. In 
this respect, our model with only one dynamical Gold- 
stone mode allows us to mark out the differences to large 
N c studies of presently used QCD low-energy models in 
a simple way. Note that the Goldstone-mode effects be- 
come even more pronounced with increasing number of 
quark flavors since the number of Goldstone modes rises 
quadratically with the number of quark flavors. 

Let us now discuss the various truncations which we 
use to study the phase transition temperature and the 
thermodynamics of our model. Our model is built around 
the standard mean-field ansatz for the effective action 
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Table I: Specification of the various truncations used to com- 
pute the equation of state. 



(large- N c ansatz), i. e. 

Z § = and Z^ = 1. 

Such an ansatz has been used extensively for studies of 
the QCD phase diagram, see e. g. Refs. [IE S3] and 
underlies also most of the recent (P)NJL studies of hot 
and dense QCD, see e. g. Refs. [H, S3 H O S3 • 

In the following we systematically extend this zeroth- 
order ansatz in two directions, namely in the number of 
derivatives and the number n of external legs of rt-point 
functions T^ n \ In order to study spontaneous symmetry 
breaking which is indicated by a non-trivial minimum of 
the order-parameter potential U(Q 2 ) in Eq. {TJ, we ex- 
pand the potential in powers of $ 2 . This results in RG 
flow equations for the mesonic rt-point functions. On the 
other hand, we perform a derivative expansion which ren- 
ders the n-point functions momentum dependent. Note 
that the latter is indispensable for a connection of the 
high- and low-momentum regime of QCD (2f| SiJ l3C| - 
see also Refs. [3lL l32j. 

Next to the zeroth-ordcr approximation (large N c ) one 
needs to include kinetic terms for the meson fields. The 
minimal truncation which allows for an inclusion of me- 
son loops is given by the so-called Local Potential Ap- 
proximation (LPA), i. e. 

Z^ = l and Z^p = 1 . 

This truncation has been used, e. g., in Refs. SH, SH, 
|35| for a study of the quark- meson model at finite temper- 
ature and density. It turns out that the LPA represents 
indeed already a major improvement with respect to the 
quality of the critical exponents. At this point, we would 
like to remark that the quality of critical exponents can 
be considered as a measure of how well the dynamics at 
the phase transition are captured. 

Next to LPA, we consider a truncation in which we 
allow for a running of the scalar wave-function renormal- 
ization, but keep the wave-function renormalization of 
the fermions fixed, i. e. 

^(fc = A)«l, d t Z^^Q, 

where t — ln(fc/A) and k denotes the IR (infrared) cutoff 
scale of our functional RG approach and A denotes the 
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UV (ultraviolet) cutoff, see Subsec. Ill Bl for details. By 
means of this approximation, vertices involving at least 
one boson field become momentum dependent. Thereby 
the quality of the results is already improved considerably 
which can be read off from the quality of the critical 
exponents 1 , see e. g. Refs. [37l. l38l. l39l] . We shall denote 
this as truncation A from now on. 

In our truncation B, we go beyond this approximation 
and allow for a running of the wave-function renormal- 
izations of the fermion fields as well, i. e. 

Z (fc = A)<l, d t Z^^0, 
Z^(fc = A) = l, dtZ^^Q. 

Such a truncation renders all vertices momentum depen- 
dent. 

At finite temperature the wave-function renormaliza- 
tions longitudinal (||) and transversal (_L) to the heat- 
bath obey in general a different RG running. In order to 
study this difference, we have to replace the part of the 
truncation ([1]) involving derivatives according to 

Z^ a i^j a + ^(d^) 2 

+z^ a i li d i ^ a + \z^ (d&f , (2) 

where we identify Z^ ^ with Z^^ in both truncation A 

and B but allow for a different running of Z^ ^ and Z^^. 
We shall denote this as truncation C. The various trunca- 
tions alongside with their labels are briefly summarized 
in Tab. [0 

As mentioned above, we expand the scalar potential 
U(Q 2 ) (order-parameter potential) in a power series of 
$ 2 up to order <I> 4 and drop all higher terms [iol [ill [42l| : 

oo r 

U(<S> 2 )^U(p,a) = J2^(p-PoT-c<r 

w A + Ai (p—po) + tt (p-paf-co , 

(3) 

where p = (l/2)<i> 2 , <&q = (oo,0) and the parameter c is 
a measure for the strength of explicit symmetry break- 
ing. Note that such an expansion represents a systematic 
expansion in n-point functions r(") where n = 2a defines 
the number of external legs and is well under control. In 
fact, the quality of an expansion of t/(<I> 2 ) in powers of 
$ 2 has been studied quantitatively in Ref. [37| at van- 
ishing temperature and for the proper-time RG in LPA 



at finite temperature in Ref. [43[ . In Ref. [ ||[ , the order- 
parameter potential U has been computed in LPA for a 
quark-meson model with two quark flavors within a RG 
approach without making use of a Taylor expansion in $. 
Note that we also solve the RG flow in the large- N c ap- 
proximation by making use of the expansion §3§ of the 
order-parameter potential U, even though it is straight- 
forward to solve the flow in this case without truncating 
the power series in Eq. ([3]), see e. g. Ref. [Hj]. In our case 
a Taylor expansion of the potential is justified since (i) we 
want to focus on the role of the momentum dependence 
of the vertices and (ii) have to ensure comparability of 
our results from our various truncations. 

From the condition dU/da = evaluated at a = a , 
we find that the RG flow of the coupling Ai and the 
minimum <Jo(k) are related: 



Al(fc)<7 (fc) = c. 



(4) 



This condition keeps the minimum at (a, n) = [ao(k), 0). 
The flow equation of the minimum ero(fc) is thus related 
to the flow of the coupling Ai in a simple way. The 
initial conditions for the various couplings in Eq. |T]) at 
the ultraviolet scale A are discussed in Sec. IIII1 



B. RG flow at finite temperature 

For our derivation of the RG flow equations of the cou- 
plings, we employ the Wetterich equation [2l|: 



1 



dtT k [ X ] - -STr 



r 



[ X ] + Rk -(d t Rk) 



with 



ri M) M = ^r fe [ x ]— 



(5) 



(6) 



where t = Ink/ A and A is the UV cutoff scale. Here, x 
represents a vector in field space and is defined as 

X T = X T (-q) ■= {yJ T (-q),Hq)^i(-q)^2(-q)) 
and 



x = x{q) ■= 



( m \ 

$>i{q) 

V Mq) J 



Thus, r£ l ' 1 ^[x] is matrix-valued in field space and so is 
the regulator function R^. In this work, we employ a 3c? 
optimized regulator function which is technically advan- 
tageous for studies at finite temperature [lj| 0, . The 
quality of such a 3c? regulator in the limit of vanishing 
temperature and chemical potential has been estimated 
by computing critical exponents of O(N) models [46j and 



1 The derivative expansion can be continued systematically b y e. g. 
including terms of the form Yj.(<I><9 M <I>) 2 in our truncation |36( . 
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comparing them to those optained with an optimized reg- 
ulator in 4d space-time [13, 5H . Details on the regular- 
ization can be found in App. [5] In the following we give 
the RG flow equations in a way which does not dependent 
on the details of our 3d regularization. Reviews on and 
introductions to the functional RG and it s appli cation to 
gauge theories can be found in Refs. [M M, \M H H3 • 
Inserting our ansatz for the effective action |T]) into 
the flow equation we find the flow equation for the 
order-parameter potential U : 



fit 17 



**_ fc 4. 



;-l/2 



1 



3 [ ^i+mj V2 

9-1/2 



-77, B (%,77v) 



1 



l + m 2 



-n F (z^,m^) >, (7) 



where w 3 = 1/(8tt 2 ), z = Z^/Z^, = z\jZ^ and 



tib{z, m) = I e j 1 / 2 t — 1 



np{z, m) = I e 5T + 1 



(8) 
(9) 



denote the Bose-Einstein and Fermi-Dirac distribution 
function, respectively. The renormalized masses of the 
bosons and fermions are given by 



1 



Z^k 2 V dp 
1 fdU 



dU n d 2 U 
2p 



dp 2 



rn„ 



Z^k 2 V dp 
h 2 p 

Tz^fk 2 ' 



2 ' #m 



(10) 

(11) 

(12) 



Here we have introduced the renormalized Yukawa cou- 
pling h = h/iZ^iZ^) 1 / 2 ). Note that the large-iV c limit 
can be deduced from Eq. ([7]) by dropping the bosonic 
loops as well as the running of the wave-function renor- 
malizations. 

In Eq. ([7]) and © we observe that the ratio of the wave- 
function renormalizations z§ = Z^/Z^ and z^ — Z^/Zjj 

effectively rescale the (dimensionless) temperature t = 
T /k in the loop diagrams. 
The anomalous dimensions, 



-ainz" 



(13) 



associated with the scalar wave-function renormaliza- 
tions longitudinal (||) and transversal (_L) to the heat- 



bath are given by 
4tj 3 



-Nr 



h 2 M<fl{t, ml) + 4 K h 4 M { fl(t, ml) 



l 4,_L 



(14) 



4wj 

z</, 



Nr 



h 2 M<fJ(t,ml) 



Anh A Mfht,ml) 



where k = Z^p^/k 2 and 
dimensions, 

% = -St In Z\ 
of the fcrmion fields read 



(15) 

Wl{Z\) 2 . The anomalous 



(16) 



2ty 



M% B l{t,m%,ml) 

+M[ F 2 B l(t,ml,ml 
Mf B ht,m%ml) 



- M l F 2 B \\(t> m l> m l S < 



(17) 



(18) 



The threshold functions M. are defined in Ref. [30] and 
App. [B] and describe the decoupling of massive and ther- 
mal modes. The RG flows of the ratios z^ — Z^/Z^ and 
z-4, = z\jZ^ read 



dtZtj, 



(19) 
(20) 



Since 77^ = 77^ and 77^ = 77^ for T/k = i = 0, which is a 
property of the threshold functions M, we have 770 = 
for t = as expected for a Poincare-invariant theory. 

The RG flow of the Yukawa coupling can be calculated 
along the lines of Ref. [3(| • We find 



d t h 2 = (4+24)h 



47Jc,/l 4 



-(FB) 



(771^,0,771^) 



-(FB),r 



y l,l 



(t,m^,, 0,m a ) 



(21) 



The threshold function L 



(FB),(4) 



1,1 



can be found in 



Ref. [30j. In the limit t = T/k — > 0, we recover 
the perturbative result d t h 2 = ((A c + 1)/(8tt 2 ))/i 4 , 
see e. g. Ref. [28]. In the chirally symmetric regime 
(m cr =m 7r ) , we find that the RG flow of the Yukawa cou- 
pling is governed by the IR attractive fixed point h* = 
for £ = T/k —>■ 0. For t — ► 00, we have /%! = <9 t /i 2 — ► 
which is due to the decoupling of the fermions for k <T. 
Note that it is special to our model with just one fermion 
species (Nf = 1) that the running of the Yukawa coupling 
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is only driven by the anomalous dimensions in the sym- 
metric regime. For example, in QCD with more than 
one quark flavor we would find additional contributions 
to the RG flow of the Yukawa coupling since the con- 
tributions from the Goldstone modes in the symmetric 
regime are not canceled completely by the ones of the 
radial mode [52|, [53|. In the regime with broken 0(2) 
symmetry in the ground state, these additional terms 
emerge when m a ^ m w . 

Let us close this section with some general remarks on 
the RG flow of ± and i at vanishing temperature. 
As mentioned above, we employ a 3d regulator function 
in order to derive the flow equations. Thus, our regu- 
larization breaks necessarily the 0{d) symmetry in the 

derivative terms of our truncation, i. e. dtZ^ , ^ d t Z^ ^ 
even for T/k — > 0. In LPA and in large- N c approxi- 
mations, this problem does not appear since the non- 
trivial momentum dependence of the propagators is ne- 
glected [H, US BUI. In studies beyond LPA, we have 
to deal with the broken Poincare-invariance at vanishing 
temperature due to the choice of our 3d regulator func- 
tion . In principle one can solve this problem by tak- 
ing care of the symmetry violating terms with the aid of 
the corresponding Ward identities. Equivalently, we can 
choose the initial conditions for the RG flow equations 
such that one finds Zjj ^ — z\ , for k — > and T — > 0. In 
the present paper, we adjust the initial conditions prop- 
erly in order to deal with this circumstance 3 . At finite 
temperature, the breaking of the 0{d) symmetry in mo- 
mentum space due to our choice of the regulator function 
does not harm our study since this symmetry is broken 
anyway. However, the choice of our regulator function 
allows us to perform the Matsubara sums analytically, 
which simplifies the numerics considerably and justifies 
our choice for the regulator function. 



III. EQUATION OF STATE AND THERMAL 
MASS 

Let us first discuss our choice for the initial condi- 
tions. Since we are interested in studying how the mo- 
mentum dependence of the propagators and vertices af- 
fects the phase-transition temperature and the thermo- 
dynamic pressure, we adjust the parameters (initial con- 
ditions at the UV scale A) of the model such that we 
obtain the same values for the IR observables for all trun- 
cations in this work. To be specific, we choose A = 2 GeV 
and /„ = = 0.030 GeV, m q = 0.088 GeV and 

Zj. j = Z\ , at k = and T = 0. Moreover we use 



Truncation 




T m ax / fn 


^max/T c 


k T =° 1 f 


large N c 


2.87 








8.96 


LPA 


2.12 


2.35 


1.11 


8.93 


Trunc. A 


1.96 


2.12 


1.08 


8.24 


Tmnc. B 


2.06 


2.26 


1.10 


8.85 


Trunc. C 


1.93 


2.12 


1.09 


8.78 



Table II: Results for the phase transition temperature T c , the 
position of the maximum in the pressure, T max , and the chiral 
symmetry breaking scale /c x sb at vanishing temperature. All 
values are given in units of the value of the order parameter 
Jtt = -y/po = (To/V2 = 30MeV at vanishing temperature. 



N c = 3 in the numerical evaluation of the RG flow equa- 
tions. Our choice for the parameters ensures that (i) 
the UV cutoff A is much larger than the physically rel- 
evant scales in order to avoid cutoff effects in the high- 
temperature regime 4 (T > T c ) and (ii) that dimensionless 
ratios of observables are close to the ones in QCD with 
Nj = 2 and Nf = 2 + 1 quark flavors. Indeed, our choice 
for the parameters yields m q / ~ 3. Since the actual 
dimensionful values for the critical temperature is of no 
interest in this paper, we give all results in units of 
This enables us to compare the results from the various 
truncations in a simple manner. 

The results for the phase-transition temperature T c for 
the various truncations are listed in Tab. |TTJ We observe 
that the phase transition temperature decreases when we 
take bosonic loops into account and resolve the momen- 
tum dependence of the vertices. Since the Goldstone bo- 
son plays a prominent role at the phase transition, it is 
intuitively clear that the inclusion of these effects affect 
the dynamics near the phase transition. Moreover, Gold- 
stone bosons tend to restore the symmetry while fermions 
tend to build up a condensate and thereby break the sym- 
metry of the ground state. However, the anti-periodic 
boundary conditions for the fermions in Euclidean time 
direction lead to a suppression of the fermionic modes in 
the vicinity and above the phase transition due to the 
absence of a zero mode. 

In Fig. [1] the influence of the bosonic fluctuations is 
visible in a different observable, namely the behavior of 
the order parameter near the phase transition which can 
be measured in terms of a universal critical exponent /?, 
/,,- ~ \T ~T C \P . We observe that the slope of the order 
parameter at the phase transition found in our large- N c 
study deviates clearly from the slope found in trunca- 
tions including bosonic fluctuations. For example, we 



2 This issue does not occur if one applies a id regulator function. 

3 From a field-theoretical point of view the adjustment of the ini- 
tial conditions of a given truncation means nothing else than 
adding appropriate counter-terms such that the theory remains 
Poincare-invariant for k — » and T — » 0. 



4 While the UV cutoff A can be chosen almost arbitrarily in LPA 
and in our large-A r c approximation, the possible choices for the 
initial conditions (UV parameters) for a given set of IR, observ- 
ables in truncations beyond LPA is restricted by the running 
of the Yukawa coupling. In this sense, these models have some 
predictive power even at vanishing temperature 
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Figure 1: (Color online) Order parameter divided by its zero- 
temperature value as a function of T/T c as obtained from the 
various different truncations for the momentum dependence. 
For all truncations beyond large- N c , the results for the order 
parameter are indistinguishable on the shown resolution. The 
different shape of the large- N c result is due to the fact that 
the value of the corresponding critical exponent differs from 
the exact value. 

find Aargc n c = 0.500 ± 0.001 and /3 LPA = 0.349 ± 0.011. 
Our large--/V c value agrees nicely with the analytically 
known value, n c = 1/2- The LPA value is in agree- 

ment with the corresponding values found in Refs. [|34| 
and [I?) • The error bars arise from a fit to our numerical 
data for the order parameter. 

Next we discuss our results for the pressure. The flow 
equation for the pressure is obtained from the negated 
right-hand side of the flow equation (J7J) for the order- 
parameter potential evaluated at the physical ground 
state $o- In Fig- LH we show our results for the pres- 
sure scaled by T 4 as a function of T/T c as obtained from 
LPA as well as from the truncations A-C; we shall dis- 
cuss the large-7V c limit below. Since we are considering 
the limit of vanishing external linear symmetry breaking 
(chiral limit, c — » in Eq. ©), we observe that the pres- 
sure approaches the Stefan-Boltzmann (SB) limit of one 
massless boson for small temperatures, p/T 4 — > 7r 2 /90 for 
T — > 0. When we increase the temperature starting from 
T = 0, we find that the pressure increases since the mass 
of the fermions, m q — hf n , decreases, see also Fig.[TJ At 
T = T c the fermions become massless while the bosons 
acquire a thermal mass m$ = m a = m^. However, the 
pressure neither takes its maximal value at T c nor does 
it agree with the Stefan-Boltzmann limit of a free gas of 
massless fermions. On the contrary, we observe that the 
pressure function reaches a maximum at T max > T c and 
approaches the Stefan-Boltzmann limit of a gas of free 
massless quarks for T > T max slowly from above. We 
would like to stress that this is a substantial difference 
to the findings in a large- N c approximations as we shall 
see below. At T — T max , we find 



P_ 



T=T„ 
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8 90 90 



Figure 2: (Color online) Pressure as a function of T/T c as 
obtained from the different truncations of the momentum 
dependence of the propagators and vertices. We observe a 
maximum in the pressure function located slightly above T c . 
For high temperatures, we observe that the Stefan-Boltzmann 
(SB) limit of a free gas of fermions (quarks) is approached 
from above. 

which is the pressure of a gas of 4A^ C massless fermions 
and 2 massless bosons. This result can be easily general- 
ized to the case with Nf fermion species. We expect that 
the maximal pressure of a QCD low-energy model with 
Nt > 2 quark flavors is given by 5 
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,N f >2 
P 



7tt2 
8 90 



4N c N f -— + N 2 f — 



Pfccc quarks 



T=T m 



„N f >2 



2Nf_ 
7 N c 



(23) 



In the second line we have divided our estimate for the 
ratio of the pressure at T = T max with the Stefan- 
Boltzmann limit of a gas of free massless quarks. 

The observation that T max > T c reflects the fact that 
chiral symmetry is not restored on all scales of the RG 
flow for T c < T < Tmax but broken on intermediate scales 
k > 0. This is depicted in Fig. [3] where we show the 
pion decay constant f„, i. e. the order parameter, as 
a function of the RG scale k in the temperature regime 
T c < T < T max . The behavior of / ff displays the in- 
terplay of the bosons and fermions in the RG flow: The 
fermions dominate the dynamics in the RG flow and tend 
to break chiral symmetry for k/T ^> 1. Below the chiral 
symmetry breaking scale fc x sB the fermions acquire then 
a mass and therefore decouple from the flow. Moreover 
all fermionic modes acquire a thermal mass ~ T. In con- 
tradistinction to the fermions, the bosonic fields have a 
thermal zero mode. Thus the RG flow is governed by the 
Goldstone boson for k/T < 1 which drives the system 



(22) 



5 Our observation of the existence of a maximum in the pressure 
function is in accordance with an estimate for the pressure func- 
tion in Ref. [33l | obtained from studying a gas of non-interacting 
quarks and mesons with temperature-dependent masses. 
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Figure 3: (Color online) RG flow of the pion decay constant in 
LPA for T/T c = 0, 1, 1.03. We observe that a regime exists 
in which the pion decay constant is finite even though the 
system is the phase with restored chiral symmetry for k — > 0. 
Note that the behavior of f n is qualitatively the same for the 
truncations A-C. 

towards the regime with restored chiral symmetry. Since 
the IR regulator scale k introduces an intrinsic length 
scale ~ 1/fc, we conclude that we have "local order" for 
T c < T < T max , i. e. the system is ordered on spa- 
tial domains 6 of length scales ~ 1/k. This implies that 
the bosons contribute still significantly to the pressure 
for these temperatures. On the other hand, the fermions 
give a significant contribution as well in this regime since 
their mass is zero for k — > 0. For T > T max , the system 
remains in the symmetric regime on all RG scales. Thus 
we have "global disorder" for these temperatures. Let us 
further discuss this phenomenon: In a Landau-Ginzburg 
type description of phase transitions, symmetry break- 
ing is indicated by a negative squared mass parameter 
(i. e. by a finite vacuum expectation value of the order- 
parameter potential). In an RG approach this parameter 
is scale-dependent. Usually we are only interested in the 
shape of the potential in the limit k — > 0. For T > T max 
the squared mass parameter of the Landau-Ginzburg po- 
tential remains positive on all scales. On the other hand, 
the squared mass parameter becomes negative on inter- 
mediate scales for T c < T < T max but becomes positive 
again in the physical limit k — > 0, see also Fig. [3] where 
the regimes of positive mass squared correspond to a van- 
ishing order parameter. Below the phase transition tem- 
perature (T < T c ) the squared mass parameter becomes 
negative at the chiral symmetry breaking scale /c x sb and 
stays negative even for k — > 0. 

Note that T c and T max differ only by roughly 10%, see 
Tab. HH Therefore the phase transition temperature can 
indeed be estimated reasonably well by just seeking for 



We would like to remark that the inverse IR regulator scale 1/k 
can be roughly associated with the side length L of the spatial 
simulation volume in lattice QCD simulations, see also Ref. [551 . 



Figure 4: (Color online) Thermal mass m$ of the scalar fields 
divided by the temperature as a function of T/T c as obtained 
from the different truncations for the momentum dependence. 
The increase of the thermal mass goes alongside with the 
decrease of the pressure function for the shown temperatures. 

the temperature for which the squared mass parameter of 
a Landau-Ginzburg-type potential stays positive on all 
scales. This strategy has been applied in Refs. [H, [13] 
in order to estimate the chiral phase transition in QCD 
from RG flows involving gluonic degrees of freedom. 

In Fig. [5]and Tab.HTl we also observe that the position 
of the maximum of the pressure is almost independent of 
the truncation. For T > T max we find that the Stefan- 
Boltzmann limit of a free gas of massless quarks is ap- 
proached faster for truncations in which the momentum 
dependence of the vertices has been resolved more accu- 
rately. In general, the behavior of the pressure function 
for T > T c reflects the behavior of the thermal bosonic 
mass m^. In Fig. 3] we show m^/T as a function of T/T c 
as obtained from LPA as well as from the truncations A- 
C. For asymptotically high temperatures we expect that 
m^/T = const, according to perturbation theory [58| . 
We find that this limit is not yet reached for the tem- 
perature range studied here. Moreover, we observe a 
difference between the results from LPA and from trun- 
cation C. The difference between the results from trun- 
cation B and truncation C, where we have also resolved 
the difference in the running of the wave-function renor- 
malizations in the Euclidean time and space directions, is 
small for T/T c < 2. This represents an a posteriori justi- 
fication for many studies where z\, = has been as- 
sumed, see e. g. Refs. (33,IH,[3l,[59j]. For a high-accuracy 
study of critical dynamics, however, our results suggest 
that higher-order derivatives may still play a role in the 
vicinity of the phase transition temperature. Therefore 
an inclusion of the full momentum dependence of the 
vertices is important for such a high-accuracy study. At 
vanishing temperature, the momentum dependence has 
been fully resolved in a study of critical exponents of 
O(N) models which has led to accurate predictions [6C| . 

Let us now discuss the results obtained from our large- 
N c approximation. In Fig. Owe compare our results for 
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Figure 5: (Color online) Pressure as a function of T/T c for 
a large-7V c approximation of our RG flow equations and for 
LPA for various ratios m^/rrig = 0, |, |. 

the pressure from a large- N c approximation in the chiral 
limit with the results from LPA for m^/mg = 0, 1/3, 1/2. 
In particular m w /m q = 1/3 resembles the situation in 
QCD with two quark flavors, where we have pions with 
a mass of roughly 130 MeV and constituent quarks with 
a mass of roughly 330 MeV. In large N c and LPA with 
m 7r/ m q > 0, we find that the pressure goes to zero for 
small temperatures since bosonic fluctuations are sup- 
pressed. With increasing temperature the pressure in- 
creases. In contrast to the results from LPA, the large- 
N c pressure function reaches its maximal value already 
at T = T c . Thus, the large- iV c results do not exhibit a 
maximum for T > T c . Comparing them with the results 
from LPA with a finite mass for the (pseudo-) Goldstone 
boson, we find that the scalar fields still contribute sig- 
nificantly to the pressure of the system for T/T c < 2, 
even for comparatively large mass ratios m x /m 9 ~ 1/2. 
This is due to the fact that the thermal contributions to 
the masses of the scalar fields dominate for T > T max . 
Therefore we conclude that in our studies the pressure 
below T c is mostly affected by the the presence of a finite 
quark mass. However, the maximum of the pressure is 
still clearly visible even though it is washed out compared 
to the chiral limit. 

Finally we comment on a possible dependence of our 
results on the UV cutoff A. We illustrate such a depen- 
dence using our results for the large-A^ approximation 
as an example, but we have checked for the other trunca- 
tions as well that our results are not spoiled by the pres- 
ence of the UV cutoff. In Fig. [HI we show the large- N c 
result for the pressure as a function of temperature for 
U/A = 0.060, 0.030, 0.020, 0.015. In the studied tem- 
perature range we find that the results for the pressure 
are strongly affected by the presence of the UV cutoff for 
fit /A > 0.020. Note that we have chosen a 3d regula- 
tor function which allows us to perform the Matsubara 
sum analytically. Nonetheless the UV cutoff still enters 
into the thermal distribution functions in the flow equa- 
tion due to their dependence on T/k. In order to 



Figure 6: (Color online) Pressure as a function of T/T c for 
various values of the UV cutoff A as obtained from the RG 
flow equations in a large- N c approximation. The parameters 
of the model for the various values of A have been adjusted 
such that the IR physics at T=0 remain unchanged. 

ensure that our results are not affected by the cutoff, we 
have chosen f„/A = 0.015 in the present work. Note that 
in the large- Af- approximation the cutoff dependence of 
the thermodynamic observables, e. g. the pressure, can 
in principle be fixed by simply adding the missing con- 
tribution (at least in the chiral limit), see e. g. Ref. [l6l |: 

6tt 2 J A e- + 1 ' 1 ' 

As soon as we allow for momentum dependent vertices 
and/or finite quark masses such a simple prescription can 
no longer be found. 

IV. CONCLUSIONS AND OUTLOOK 

We have discussed the truncation dependence of the 
thermodynamics of a simple model consisting of strongly 
interacting relativistic bosons and chiral fermions. Our 
results indicate that the derivative expansion is a good 
and valid approximation for a study of the thermal prop- 
erties of Yukawa- type theories, such as the pressure or 
the behavior of the order parameter below the critical 
temperature. We find that the pressure exhibits a max- 
imum above the critical temperature when we include 
bosonic fluctuations. The origin of the maximum is re- 
lated to the fact that the thermal masses of the scalar 
fields are still small in a small window above the criti- 
cal temperature and therefore contribute significantly to 
the equation of state. For very high temperatures, the 
pressure then approaches the Stefan-Boltzmann limit of 
a gas of free massless fermions from above. The conver- 
gence to the Stefan-Boltzmann limit is faster the better 
the momentum dependence of the vertices is resolved. 

On the other hand, we have also shown that the critical 
temperature shows a stronger dependence on a variation 
of the truncation. In order to clarify whether our results 
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for T c can be considered as having almost converged, we 
eventually need to fully resolve the momentum depen- 
dence of the propagators and vertices as it has been suc- 
cessfully done in Refs. [H HU, M, M, H, H M, at 
vanishing temperature. This is deferred to future work. 

Whether the present findings, in particular the behav- 
ior of the pressure above T c , are visible in the thermo- 
dynamic pressure of QCD can of course not be answered 
by the model used in this study. First of all, the model 
is not confining and therefore the pressure below T c is 
overestimated. Within a functional RG approach, this 
can be resolved by combining the findings of recent func- 
tional RG studies of QCD including gauge degrees of free- 
dom [H, [3(J [H, ■ This would allow us to compute the 
pressure from first principles in the sense that it would 
rely on only a single input parameter, namely the value 
of the strong coupling at the initial RG scale. Above the 
critical temperature, the pressure is mainly dominated by 
the presence of gluons. Still, the contribution from the 
mesons in the vicinity of the phase transition can be es- 
timated from Eq. (|23|) to yield a 10% correction in QCD 
with two massless quark flavors. Therefore it might be 
possible to observe this effect in future lattice or func- 
tional RG studies of the thermodynamic pressure. 
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Appendix A: DECOMPOSITION OF 
MOMENTUM TENSORS AT FINITE 
TEMPERATURE 

Throughout this paper we encounter integrals of the 
following type 



defined as follows: 



Z(T) = T J2 



jd-l 



(2n) 



(Al) 
(A2) 



where T denotes the temperature, po is a discrete variable 
depending on an integer n and the function /(.,.) sim- 
ply represents a placeholder for a general function which 
is symmetric under p^ — > — p^. Due to symmetry, the 
tensor structure can be simplified as follows: 



PvPv -> Pl P \v 



P 2 



P ± 
d-1 » v 



(A3) 



where n„ = (1,0) denotes the heat-bath velocity and the 
transversal (_L) and the longitudinal (||) projectors are 



P. 



and 



p, 



fiiy 



(A4) 



Appendix B: THRESHOLD FUNCTIONS 

The regulator dependence of the flow equations is con- 
trolled by (dimensionless) threshold functions which arise 
from Feynman graphs, incorporating fermionic and/or 
bosonic fields. Let us first introduce the so-called di- 
mensionless regulator-shape function r^(x) and r$(x) for 
bosonic and fermionic fields. These functions arc implic- 
itly defined by the regulator function Ri as follows: 



R*(po, {Pi}) = Zj(P0, {Pr})p 2 r B {p 2 /k 2 ) 
for the bosonic fields and 

RipiPoAPi}) = z ^,k(Po,{Pi})P r 4>(p 2 /k 2 ) 



(Bl) 



(B2) 



for the fermionic fields. 

In this work, we employ a 3c? optimized regulator-shape 
function [HIM SI : 



and 



r^(x) 



r^(x) 



- - 1 ) 6(1- x) 

x 



1 0(1 



(B3) 



(B4) 



In the following, we use these regulator shape functions 
whenever we evaluate the integrals and sums in our gen- 
eral definitions of the threshold functions. Those thresh- 
old functions which are not defined below can be found 
in Ref. [H. 

In order to define the threshold functions, it is conve- 
nient to define dimensionless propagators for the scalar 
fields ((f)) and the fermions (ip), respectively: 



G^(xl,u}) 



1 



(B5) 



where x, 



o 



2-Kni and x^ 



z^Xq + x(l + r$) +uj 

and G^(xt,u) = — t — , (B6) 

z^Xq + x(l + r^y + oj 

2n(n + 1/2)1 Note that the 
Matsubara modes are dressed with the ratio of the wave- 
function renormalizations longitudinal and transversal to 
the heat-bath, z$ = Z\/Z& and z^ = z\/Z^. Moreover 
it is convenient to define the dimensionless propagators as 
obtained from the application of the optimized regulator- 
shape function defined in Eqs. (|B3[) and (IB4|) : 



1 



(B7) 



1 + Z^Xq + UJ 

and G^x^uj) = — t . (B8) 

1 + z\xl + UJ 
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For the wave-function renormalization of the scalar 
fields, we find: 



°° f°° d-3 ~ { 

= (d-l)t ^2 / dxx~d t < x{l+r^)G^{x^,uj^)x 



n——oo 
72 



d+ 1 / d 



44^(4)^) 



2x ( d n i i> \ 
+ (-^G^(x^,u>,p) 



(B9) 



Employing the optimized 3d regulator function specified 
above, we obtain 



M { 4 F l(t,u;) 

= 4* J2 (i + 44)(^(4^v)) 

n— — 00 ^ 

^ (&/>(4>«v)) - (^z^+J) (&/.(4> w v>)) 



+ 



3g (^(^0'^)) + (^(^O 



d-3 



(BIO) 



The threshold function for the corresponding longitudi- 
nal contribution reads: 



00 f°° d 3 f 
t ^ y dxx~9 t < a;(l + r^,)G^(xQ ,w^,)x 



2x 



1 (rf4) 2 



(1 + r,p)G^(xg , uty) 



+ — — (! + f>)GV(4,<*ty) 
\ da; n 



+ 44£v>(4> w V>) 



{dxtY 



+3 



dXg 



n— — 00 \ / V. 

+(26^4 + 6) (&/,(4,c^)) 4 

]4(i + 4^0) (^(4, ^)) 5 } • (bii) 



-325.1 



Apart from these contributions, the scalar wave-function 
renormalization also gets a contribution from a fcrmion 
loop which is proportional to the vacuum expectation 
value of the scalar field: 

I f°° d 3 

M^(t,w^,) = -z^t ^2 J dxx^d t G^{x^,uj^) x 



d 2 



+2 44 L dx ^ 2 ^( x o.^) 



1 'i-^u* 



\ / n— — 00 

x {3 - 164a#S*(*o.aty)} . (B12) 

Finally, there is a contribution to the scalar wave- 
function renormalization from a pion-sigma loop which 
is proportional to the four-boson coupling: 

= ~8~ E / dxx^d t G<j>{xt-,uJ4,,i) (Gir(x$,u) a )^ x 

X 18^4^(4'^) _ 2 } + ^ ^) 

x (g (4, W(T )) 2 [2^(4,^+^(4,^) 
-4^4^0(4' w tt)(^(4' ^o-)) [^(4i ^) 

+3 g (4, u> a j\ I + K ^ . (B13) 



For the fermionic wave- function renormalization longitu- 
dinal to the heat-bath, we find 



= 2 E I dxx^d tJI ^- ) l-G^xt,^)x 



x T^^(4^0) 

\ dXn 
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d-1 



. ^ \ ( r i ^ \ \ 2 tions can be performed analytically. Since the resulting 

x o yip\ x o ' W V>J \y<t>\ x o ' W 4>JJ x expressions are rather complicated and of no interest for 

the discussions in the present paper, we do not show them 

I G. . i . i r 2 ! I — ! G. .(.n. ■ u.-. .) 1 . 



^ | n ip o / i ^ i r V> \l here. However, for the numerical evaluation of the RG 

1-— \g^{x Q ,uj^) + 2\ l- — j ^(Zq , U4,) j> . flowg wg haye made uge of thg fact that the Matsubara 

sum in these functions can be performed analytically. 



(B14) 

Note that the Matsubara sums in the threshold func- 
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